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ABSTRACT 

Several luminosity indicators have been found for Gamma- Ray Bursts (GRBs) 
wherein measurable light curve and spectral properties are well-correlated with 
the peak luminosity. Several papers have each applied one different luminosity 
relation to find redshifts for BATSE GRBs and claim to identify specific bursts 
with z > 8. The existence of such high redshift events is not surprising, as BATSE 
has enough sensitivity to see them and GRBs are expected out to the redshift 
of the first star formation. To improve results we used five luminosity relations 
with updated calibrations to determine redshifts with error bars. Combining 
these relations, we calculated the redshifts of 36 BATSE GRBs with claimed 
z > 8. Our results include 13 bursts with our derived best redshift Zb es t > 8, 
which looks promising at first. But the calculated redshift uncertainties are 
significantly large in these selected cases. With only one exception, all of our 
bursts have z\ a [ ow < 9. The one exception (BATSE trigger 2035) is likely a short 
duration burst at z > 4. Our best case for a very high redshift event is BATSE 
trigger 3142 with Zb es t > 20 and z la i ow = 8.9, however we can only say z > 4.1 at 
the two-sigma confidence level. In all, we cannot point toward any one BATSE 
burst as confidently having z > 8. One implication is to greatly weaken prior 
claims that GRBs have a steeply rising rate-density out to high redshifts. 

Subject headings: Gamma-Ray: Bursts - Early Universe 



1. Introduction 

Long duration Gamma-Ray Bursts (GRBs) are expected to be visible to high redshifts. 
This is because they are caused by the core collapses of massive stars, which have very short 
lifetimes (Woosley & Bloom 2006), so the first bursts in our Universe should be shortly after 
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the first star formation. The WMAP detection of reionization at redshifts 20+g° (Kogut 
et al. 2003), 17 ± 4 (Spergel et al. 2003), or ~ 10 - 20 (Spergel et al. 2006) set the 
scale for the redshift at which the first significant star formation occurs. On the reasonable 
assumption that the GRB rate is proportional to the massive star formation rate, models of 
star formation in our Universe can be translated into predictions as to the relative rate of 
GRBs as a function of redshift. Detailed calculations show that the first GRBs should have 
redshifts of > 20 (Bromm & Loeb 2002; 2006) or ~ 60 (Naoz & Bromberg 2007). GRBs 
have already been identified at z — 6.295 (Totani et al. 2006) and z = 6.6 (Grazian et al. 
2006; Tanvir et al. 2006). The lack of bursts with higher spectroscopic redshifts is likely due 
to the severe observational difficulties in recognizing the events and measuring their spectra. 
Detailed models for the fraction of Swift bursts with z > 8 give 10% (Bromm & Loeb 2002), 
3% (Le & Dermer 2007), and ~5% (Bromm & Loeb 2006). 

The BATSE detectors onboard the Compton Gamma Ray Observatory have a deep 
sensitivity threshold and most of the sky was monitored for nine years, so there are likely to 
be many high redshift GRBs hidden in its catalog. BATSE can detect typical GRBs out to 
redshifts of 20-30 (Lamb & Reichart 2000). The fraction of BATSE bursts with z > 8 should 
be 4% (Bromm & Loeb 2002). With 1637 bursts in the fourth BATSE catalog (Paciesas et 
al. 1999), we expect ~ 70 z > 8 GRBs. 

One way to recognize high redshift BATSE bursts is to use the various luminosity in- 
dicators. Just as with Cepheids and supernovae, readily observable light curve and spectral 
properties are well-correlated with the luminosity, so that the observed brightness plus the 
inverse- square law will allow a luminosity distance to be derived. For some fiducial cosmology 
with a distance/redshift relation, the redshift can then be calculated. The luminosity indi- 
cators include the spectral lag time (Norris, Marani, & Bonnell 2000; Norris 2002; Schaefer 
2004), the variability (spikiness) in the light curve (Fenimore & Ramirez- Ruiz 2000; Reichart 
et al. 2001; Meszaros et al. 2002; Kobayashi, Ryde, & MacFadyen 2002; Li & Paczynski 
2006), the photon energy of the spectral peak E pea k (Schaefer 2003; 2007; Yonetoku et al. 
2004), the minimum rise time in the light curve (Schaefer 2007), and the number of peaks 
in the light curve (Schaefer 2007). A very-high-luminosity burst is recognizable as having a 
near-zero lag, a very spiky light curve, a high E peak , a very fast rise time, and many peaks. 
The luminosity relations (i.e., luminosity as a function of the luminosity indicator) are cal- 
ibrated with 69 GRBs that have spectroscopic redshifts (Schaefer 2007). The luminosity 
relations vary substantially in their accuracy, and one of them (the number of peaks) only 
provides a lower limit on the luminosity. 

Various previous researchers have applied one of the luminosity indicators to a large 
number of BATSE bursts (Fenimore & Ramirez- Ruiz 2000; Schaefer, Deng, & Band 2001; 
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Lloyd-Ronning, Fryer,, & Ramirez- Ruiz 2002; Band, Norris, & Bonnell 2004; Yonetoku et 
al. 2004; Kocevski & Liang 2006). The goal was to measure the redshifts of BATSE bursts 
for various demographic reasons. One such purpose was to measure the burst event rate as 
a function of redshift so as to derive the evolution of the universal star formation rate, while 
another purpose was to determine the burst luminosity function and its evolution. These 
studies all identified individual bursts that apparently had very high redshifts. 

Band, Norris, & Bonnell (2004) used the lag-luminosity relation to derive the redshifts 
to 1194 BATSE bursts, for which they identified 44 bursts with z > 8 (4%) and 17 bursts 
with z > 12 (1.4%). Yonetoku et al. (2004) used the -E pea fc-luminosity relation to derive 
the redshifts of 745 BATSE bursts, for which 77 have z > 8 (10%) while none have z > 
12. Fenimore & Ramirez-Ruiz (2000) used the variability-luminosity relation to derive the 
redshifts of 220 BATSE bursts, finding 38 bursts with z > 7 (17%) and 19 bursts with 
z > 12 (8.6%). If correct, these very high redshifts would be of high interest and useful for 
cosmology. 

Unfortunately, as these researchers were well aware at the time, the highest derived 
redshifts have substantial uncertainties. To give an example, consider the lag-luminosity 
relation, where the burst luminosity is inversely proportional to the lag. The highest lumi- 
nosity events will all have small lags, but the measurement uncertainty of lags is roughly a 
constant (depending on the number of detected photons and the light curve shape) so it will 
be impossible to distinguish between a very high luminosity burst (with a very short lag) at 
a very high redshift and a moderately high luminosity burst (with a moderately short lag) 
at a moderately high redshift. To illustrate this, consider BATSE trigger 2850 for which the 
lag was measured to be 0.008 seconds, but for which the one-sigma range is between -0.044 
and 0.092 seconds (Band, Norris, Bonnell 2004). The best fit value for the lag happened 
to be near zero out of a large uncertainty range and this was followed through to yield a 
formal best fit redshift of 22.1. But the one-sigma range of lag runs to 0.092 sec, and this 
corresponds to a luminosity 11.5 times fainter and a redshift of ~ 0.1. Thus, in the ordinary 
course of Gaussian errors, many moderate redshift bursts with a lag of, say, 0.1 sec will 
appear at very high redshift due to measured lags of 0.01 sec or even 0.001 sec. As such, 
normal measurement errors will smear the burst redshift distribution so as to artificially 
populate the highest redshifts. Thus, keeping track of the uncertainties is important. Only 
the Yonetoku et al. (2004) paper actually quoted their redshift uncertainties. Of their 77 
bursts with the best estimate of the redshift as z > 8, only 9 have the one-sigma range above 
z = 8, while only one has its two-sigma range barely requiring z > 8. 

Another challenge for the idea of these BATSE bursts being at very high redshift is 
that derived redshifts from multiple luminosity indicators should agree that the redshift is 
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high. For example, it would not be convincing to have a redshift of 12 ± 1 derived from 
the burst's lag if the variability and E peak values are pointing to z ~ 3. However, if all the 
luminosity indicators unanimously pointed to z ~ 12, then we would have some confidence 
that the redshift is indeed very high. To date, no such cross-comparison has been reported. 
The use of multiple indicators has a further advantage that the independent measures of the 
luminosity can be combined as a weighted average so as to substantially reduce the size of 
the derived one-sigma region. 

Another problem with the old reported BATSE redshifts is that they were made with 
old definitions and calibrations. Many more GRBs with spectroscopic redshifts are now 
available for calibrating the luminosity relations, and realistic systematic errors are now 
known (Schaefer 2007). 

This paper reports on our concerted attempt to identify very high redshift GRBs in 
the BATSE database. In particular, we tested whether any of the specific bursts reported 
to have very high redshifts can indeed be confidently concluded to be at very high redshift. 
For this, we are using five luminosity indicators so as to look for consistency. We are using 
the latest definitions and calibrations for the luminosity relations. We are keeping track of 
the uncertainties in each indicator so as to get a best fit redshift as well as a quantified 
one-sigma range. The result is that we can confidently test whether the previously identified 
very-high-redshift BATSE events are really at such high redshifts. 



2. GRB Selection and Properties 

Initially, we looked through the published lists of BATSE burst redshifts for cases where 
one burst was consistently quoted to be at very high redshifts. No such cases were found. 
This in itself is a foreshadowing of our final answer. Nevertheless, we soon found a variety 
of the usual problems concerning consistencies between measurements and old calibrations. 
Thus, we realized that we would have to derive many of the light curve properties on our 
own and then apply our recent calibrations. 

We decided to look at the 36 BATSE bursts with the highest chances of being at 
very high redshifts as based on prior publications (Band et al. 2004; Yonetoku et al. 2004; 
Fenimore & Ramirez- Ruiz 2000). These three prior publications quote independent redshifts 
for hundreds of BATSE bursts based on one of three GRB luminosity indicators (lag, E peak , 
and variability respectively). We selected the 12 GRBs with the highest redshifts from each 
of the three sources. (Again, the lack of overlap in these selections point to a lack of confident 
high-z bursts.) Two bursts originally chosen did not have any available E peak , so the next 
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two highest redshift bursts from the same source were selected in their place. The 36 selected 
bursts all have z > 8 as reported by one of the three papers. This selection criterion will 
inevitably lead to bursts with one extreme indicator, discrepancies between indicators, and 
large error bars in redshift. 

Table 1 contains the measurable burst properties leading up to the derived bolometric 
peak flux. Column 1 contains the BATSE trigger number that identifies each burst. A 
footnote identifies which of the papers was used for selecting the burst. Column 2 lists the 
T 90 of each burst, with this duration measuring the interval in seconds between the time 
when 5% to 95% of the total counts have been detected. Column 3 has the observed peak 
flux over a one second time interval in the energy range of 50-300 keV (P50-300 with units of 
photon cm -2 s _1 ). Both P50-300 an d T 90 have been taken from the Fourth BATSE catalog 
and the 'current catalogf] (Paciesas et al. 1999). GRB spectra can be reasonably described 
as smoothly broken power laws, as given by the Band model (Band et al. 1993). In this 
model, the low energy power law index is a, the high energy power law index is /3, and the 
photon energy of peak spectral flux is E pea k. Columns 3-5 contain the values for E pea k, a, 
and (3. These values are those reported by Band, Norris, & Bonnell (2004)@. Of these, the 
majority are from unpublished measures by Robert Mallozzi (c.f. Mallozzi 1995). The values 
placed in square brackets are typical default values adopted when the value itself was not 
measured. The spectral parameters and the peak fluxes were used to find the bolometric 
peak fluxes (Pboio) as in equation 6 of Schaefer (2007). The value was calculated for the best 
redshift (see Table 4), although the dependancy on the assumed redshift is relatively small. 
The resulting Pboio (i n column 7) is over a one second time interval and has units of erg cm -2 
s _1 . The uncertainty of Pboio is found by propagating the uncertainties of P 50 _3oo, Epeak, 
and j3. 

Table 2 lists the luminosity indicators for each burst. {E pea k is also a luminosity indi- 
cator, although it is tabulated in Table 1.) Column 2 contains the lag time (ri ag ) of a GRB 
which is the time in seconds between the hard (100-300 keV ) and soft photon (25-50 keV) 
light curves (Norris, Marani, &; Bonnell 2000). The values listed are taken from Band et al. 
(2004), with the quoted error bars being for the positive lag direction. These values agree 
within error bars with the lags reported for these same bursts by Hakkila et al. (2007). If 
Band et al. (2004) found a negative or zero lag, we took the lag to be 0.001 seconds. Column 
3 has the variability (V) which is dimensionless and measures the spikiness of the burst light 
curve. It is measured by calculating the variance of the light curve from a smooth version 



1 http://cossc. gsfc.nasa.gov/docs/cgro/batse/BATSE_Ctlg/index. html 



2 See also 



http://cossc.gsfc.nasa.gov/docs/cgro/analysis/lags/ 
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the light curve (Fenimore & Ramirez-Ruiz 2000). We used the BATSE 25-300 keV light 
curves with the definition of variability as given in equation 9 of Schaefer (2007). Column 4 
contains the minimum rise time {ruTmin) which is the shortest amount of time in seconds for 
the light curve to rise by half its peak. The number of significant peaks (N pea k) in each burst 
light curve is listed in column 5. A peak is counted if its peak flux lies above the background 
is at least 25% of the largest peaks flux and there is a dip between it and the surrounding 
peaks of at least 25% of the peak flux. We have calculated V, r RTmin , and N peak from 0.064 
ms concatenated BATSE light curves for 25-300 keV@. 

3. Luminosity Relations 

The luminosity relations are calibrated equations that give the burst's peak luminosity 
(L) as a function of the measured luminosity indicators. To date, eight different luminosity 
relations have been proposed and verified. In this paper, we will be using only five of these 
relations, the same five as used by Schaefer (2007) for constructing the 69 GRB Hubble 
diagram. The E pea k — E 1 relation of Ghirlanda, Ghisellini, & Lazzati (2004) was not used 
as the required jet breaks in the afterglow light curves were not seen for BATSE bursts. 
The E pea k — E ison relation of Amati (Amati et al. 2002; Amati 2006) was not used as its 
luminosity indicator is duplicated in the E pea k — L relation which we are already using. The 
E P eak — T 0A5 — L relation of Firmani et al. (2006) is not used as its input and physics is 
similar to that of the E pea k — L relation which we are already using. 

The five luminosity relations that we use are summarized in Table 3. The lag-luminosity 
relation is fairly tight, but the bunching at near-zero lags makes for difficulties in distinguish- 
ing high-luminosity and very-high-luminosity bursts. The variability-luminosity relation is 
one of the poorest of the eight in terms of accuracy, but there is a substantial spread in 
variability between the high-luminosity and the very-high-luminosity bursts. The E pea k- 
luminosity relation was originally predicted theoretically and then observationally confirmed 
(Schaefer 2003), with later verification (Yonetoku et al. 2004). The r fiTmin -luminosity rela- 
tion was also predicted and confirmed in 2002, although this relation is also relatively poor 
in accuracy. The iV pea fc-luminosity relation is another prediction from 2002, and even though 
it only provides a limit on the luminosity, the relation works in the right direction for the 
purposes of this paper as a many-peaked bursts must have very high luminosity. 

The luminosity relations must be calibrated with bursts of known redshift (and hence 
with known peak luminosity). Schaefer (2007) has collected 69 GRBs for use in calibrating 



3 Available at |http://cossc.gsfc.nasa.gov/docs/cgro/batse/batseburst/sixtyfour_ms/index.litml 
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the five luminosity relations, and we have adopted these calibrations (see Table 3). Impo- 
rantly, these calibrations are for the same definitions of luminosity indicators as we use in this 
paper. So, for example, we should not use the V values given by Fenimore & Ramirez-Ruiz 
in the equation quoted in Table 3. 

A possible concern is that the luminosity relations are calibrated with a variety of 
different satellites and then applied to BATSE. For light curve luminosity indicators, the 
only systematic difference could arise from the use of different energy bands. But here, 
the chosen bands are closely similar (see Table 5 of Schaefer 2007), so systematic problems 
should be negligible. For the E pea k relations, the different satellites are all measuring the 
same quantity without any known biases, so as long as the (possibly asymmetric) error bars 
are propagated, the joint calibrations should apply to any individual satellites data. The 
real proof that the calibrations are not satellite dependent is that Schaefer (2007) found the 
best fit calibrations for each satellite to be identical to within the error bars. 

The accuracy of the derived redshifts from combining all five luminosity relations is 
substantially better than can be obtained by any one relation alone (Schaefer 2007). The 
overall accuracy can be quantified by comparing derived redshifts versus the spectroscopic 
redshifts. For 69 GRBs, the one-sigma scatter about the known redshifts is 26% (Schaefer 
2007). And the reduced chi-square of this comparison is close to unity, with the implication 
that the calculated error bars are realistic. This is the basis for knowing that the luminosity 
relations are reliable and accurate to within the stated error bars. 

The luminosity relations do not apply to short duration bursts, presumably because they 
have a different physical mechanism. The dividing line between long-duration and short- 
duration bursts is around 2 seconds as measured in the BATSE rest frame (Kouveliotou et 
al. 1993). Presumably, such a division is better done in the rest frame of the burst, which 
equals 2/(1 + z) s, where z is the typical redshift of BATSE bursts. For the typical redshifts 
of BATSE bursts of roughly 1, the division between long and short should be roughly 1 
second in the burst rest frame. 



4. Redshifts 

Each of the five luminosity relations will yield a luminosity with its error bar. For each 

i 

of these luminosities (along with Pboio) we can derive a luminosity distance as (L /[47r Pboio}) 2 ■ 
Then, we can get a distance modulus and its uncertainty (fi ± a^) for each of the 36 GRBs. 
But each of the luminosity indicators must be corrected to the rest frame of the burst. So 
we have set this up by evaluating jj, and cr^ for all redshifts (each corresponding to a distance 
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modulus fj, z ) from 0.01 to 20 at an interval of 0.01. For each luminosity relation and each 
trial redshift, we can evaluate the chi-square as [(// — ji z )/a^)} 2 . The resultant chi-squares 
for all five luminosity relations are then coadded to get a total chi-square which is a function 
of redshift. Our best redshift (zb est ) is the value of z for which the chi-square is the smallest. 
Our one-sigma confidence region for the derived redshift is the range of redshift over which 
the chi-square is within 1.0 of the minimum chi-square. Our two-sigma region is where 
the chi-square is within 4.0 of the minimum. For each individual luminosity relation, we 
can similarly derive a best redshift (zi ag , zy, ZEpeak, zrt, and z^peak) and the one-sigma 
confidence interval. These are listed in Table 4. The error bars are generally asymmetric, 
so we have represented the one-sigma confidence interval in parentheses immediately after 
the optimal value. Some of these redshifts or their confidence ranges are open-ended to high 
redshift, and we represent this by capping the upper limit at a redshift of 20. 

How consistent are the five derived redshifts from the five luminosity relations? A 
quantitative way to answer this question is to look at the chi-squares for Zb es t- With three 
degrees of freedom for our chi-square fit (the limit from N pea k is never constraining), we 
expect that the average chi-square will be around 3, while a three-sigma discrepancy would 
require the chi-square to be greater than 3+9=12. For our 36 bursts, we find that the typical 
chi-square is around 3 and the largest chi-square values (for BATSE triggers 2345 and 7390) 
are near 7. From this, we conclude that all of our luminosity indicators are giving consistent 
redshifts. 

The one-sigma uncertainty ranges on Zf, est can be fairly large. One reason for this is 
that the very-high luminosity case gives near-zero indicators for both lag and risetime, so 
that the usual measurement uncertainties only put a lower limit on the luminosity. Also, the 
high E peak values are difficult to measure with accuracy. Another reason is that the slope of 
the redshift versus luminosity distance relation changes at high redshifts, so that the upward 
error bars in redshift extend over a larger numerical range. As such, our expected large 
redshift uncertainties for the 36 bursts in this paper are not characteristic of normal bursts. 

In 13 cases, the z^ est value is > 8, which would be of interest, but the one-sigma range 
extends down to moderate redshift. For example, BATSE trigger 1734 has a best redshift 
of 14.5 but a one-sigma confidence region of > 8.3 and a two-sigma confidence region of 
> 4.1. In such a case, we cannot confidently identify the burst as z > 8, because the redshift 
could be smaller at only the 1.1-sigma confidence level. That is, the redshift of trigger 1734 
can easily be 4- 6, which is no longer novel. For convenience, we have added a last column 
to Table 4 which lists the redshifts of the lower edges of the one-sigma confidence regions 
(ziaimv)- Let us examine the redshifts for four cases so as to see the range of our results. 

The first case (BATSE trigger 111) was originally selected because its E pea k value was 
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sufficiently high for Yonetoku et al. (2004) to quote that z Epe ak = 11.2, whereas we get a 
more moderate ZE pea k — 5.5. This difference is not due to the E pea ^ value, as that from 
Band et al. (2004) is close to that used by Yonetoku et al. Rather, the difference is mainly 
caused by the difference in the slope of the E pea k — L relation, where they used 2.0 ± 0.2 
while we used 1.68 ± 0.05. (The two slopes for the relation are different because Yonetoku 
et al. [2004] only used 11 GRBs for calibration, whereas Schaefer [2007] used 64 GRBs for 
calibration.) In any case, Yonetoku et al. gave the lower edge of their one-sigma range for 
ZEpeak as 4.72, so there is no significant discrepancy. In the meantime, all the other four 
luminosity relations point to a moderately low redshift, as the burst has a long lag, smooth 
light curve, slow rise time, and only one peak. The Zb es t is 1-4, which is only just outside the 
one-sigma range for ZEpeak- In this case, ZE pea k has a large deviation from zt, es t, and this is 
not surprising because the burst was selected for its extreme ZE pea k- 

The second case (BATSE trigger 493) was originally selected because Yonetoku et al. 
reported that ZE pea k = 11.9, even though the lower edge of their one-sigma range was to 4.6. 
We get ZE P eak = 4.3 with a somewhat smaller value for E pea k. All five individual redshifts 
are consistent with Zbest = 3.7. For this redshift, the duration in the burst rest frame will 
be T 90 /(l + z) — 1.0 s, which still allows the event to be a long duration GRB. Our final 
redshift has a one-sigma range of 2.3-6.1, so we are confident that this GRB is not at high 
redshift. 

The third example is our best case for a very high redshift GRB, BATSE trigger 3142. 
Our Zbest is > 20, which appears to be of interest. However, Z\ a i ow = 8.9 while a two-sigma 
limit is > 4.1, so a reasonable idea is that the burst has an 'unexciting' redshift. Nevertheless, 
the burst has a very short lag, a high variability, a high E pea k, a very fast rise time, and 8 
peaks, so this remains a good candidate for a very high redshift burst. 

The fourth case is BATSE trigger 2035, which has a z^est > 20 and Z\ a i ow = 14.9. This 
burst has a very short lag, a very high variability, and the shortest of the rise times, all 
of which point to a very high luminosity GRB. When combined with the faint Pboio, these 
three indicators point to z > 10. However, BATSE trigger 2035 is a burst whose duration 
is near the low end of the range for long-duration bursts, with T 90 = 5.1 s. If the burst were 
indeed at z — 14.9, then the time dilation of the light curve would imply that the duration 
was 0.3 s in the burst's rest frame and that the event would then be a short- duration burst. 
But the luminosity relations do not apply to short bursts, so we could not then derive the 
high redshift. Thus, we have a logical inconsistency with BATSE trigger 2035 being a long- 
duration GRB with z > 4. Our best solution is to say that the event is actually a short 
duration burst with z > 4 or so. In any case, this event cannot be a very-high-redshift 
long- duration GRB. 
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If BATSE trigger 2035 is a short duration event, then its redshift must be z > 4 so that 
its T 90 duration would have been less than one second in its rest frame. Is it plausible to have 
a short burst at such high redshifts? The reason for asking this is that the first discoveries 
of optical afterglows of short bursts were associated with galaxies at low redshifts (z = 0.16 
for GRB050709, Fox et al. 2005; z = 0.258 for GRB050724, Prochaska et al. 2005; and 
z = 0.547 for GRB051221, Soderberg et al. 2006), and this suggested that all short GRBs 
are at low redshifts. However, the lack of any host galaxies to very deep limits in five short 
GRB error boxes (Schaefer 2006) indicates that many of these events are at moderate or 
high redshifts. Berger et al. (2007) later obtained similar limits and concluded that many 
short bursts are at z > 1. Also, from a theoretical population synthesis approach for double 
neutron star binary mergers, a high proportion of short GRBs are predicted to occur at 
z > 1 (O'Shaughnessy, Belczynski, & Kalogera 2007). In all, it is reasonable to think that 
BATSE trigger 2035 is a short duration GRB with z > 4. 

This same duration dilemma applies to 7 of our 36 bursts. That is, the bursts with 
^90 ^ 9 s cannot be identified as z > 8 events because then their durations in their rest 
frames would be under 1 second and they would not be long duration GRBs for which the 
luminosity indicators apply. A possible solution is that these are long duration GRBs with 
moderate redshift or short duration bursts with moderate redshift. Only 3 of our bursts 
(BATSE triggers 2035, 2850, and 2944) have T 90 /(l + z laiow ) < 1 s. These events could 
be short duration bursts at moderate redshifts. Or they could be long duration bursts 
whose moderate redshifts are only somewhat smaller than z lcr / ow . This last possibility is not 
surprising, as we expect 16% of the bursts to have their real redshift below z 1(T i ow , and the 
fraction might be substantially higher due to our selection of bursts for extreme- ,2 by only 
one indicator. 

The goal of our research was to test the claims to identify specific BATSE bursts with 
very high redshifts. When we began this work, we were hoping to be able to identify GRBs 
with z > 8. Indeed, 13 out of our 36 GRBs have Zb est > 8. However, our candidate 
very-high-redshift bursts all have too large of error bars to decisively say they are at high 
redshifts. With one exception (BATSE trigger 2035, see above), all bursts had Z\ a \ ow < 9, 
so we cannot confidently state that any of them are of very high redshift. Some of these 
candidates might be at z > 8, or alternatively they all might be at 'unexciting' redshifts 
of < 6. The less interesting alternative might even be preferred because moderate-redshift 
events are apparently more common than high-redshift events so that normal observational 
scatter in the derived redshifts will have the high-z best population dominated by moderate- z 
bursts. 



5. Discussion 



We have not been able to find any GRBs that are confidently at z > 8. But we know 
that roughly 4% of the BATSE catalog should be at comparably high redshifts (Bromm & 
Loeb 2002). So where are they? Part of the answer might be simply that our Zb es t > 8 bursts 
are indeed at z > 8, although we cannot prove that because the uncertainties are too large. 
The other part of the answer is likely simply that the z > 8 GRBs are in the fainter half 
of the BATSE catalog. These fainter bursts are too faint for Band et al. to have derived 
their lags, too faint for Fenimore & Ramirez-Ruiz to derive their variabilities, and too faint 
for Yonetoku et al. to measure their E pea ^. As such, the faint half of the BATSE catalog 
has not even been examined for high redshifts. It is not clear whether the majority of the 
very high redshift events should appear in the faint half of the BATSE catalog (as befits 
their greater distance) or should appear in the brighter half of the BATSE catalog (as befits 
the evolution of the GRB luminosity function and the greater volume at high redshift, of. 
Schaefer 2000). In any case, any very high redshift GRBs in the faint half of the BATSE 
catalog will have only larger uncertainties than those in the brighter half. The solution as 
to where are the high redshift BATSE bursts is not likely to depend on the various types 
of Malmquist biases, as Schaefer (2007) has made detailed calculations on the Malmquist 
and gravitational lensing biases only to find that all the biases together nearly cancel out 
resulting in negligible effect. In all, we conclude that BATSE catalog could well include 
many z > 8 GRBs, but that the redshift uncertainties will always be too large to allow for 
the confident identification of any such GRB. 

What does all this say about the many demographic studies which extend to z > 8? 
The primary problem here is that the normal scatter of derived redshifts will appear to move 
many moderate- and high- redshift GRBs out into the very-high-redshift regime. The large 
number of lower- z bursts will produce many more Zbest > 8 events than there are z > 8 
bursts that are made to appear of low or moderate redshift. Thus, the number of apparent 
very-high- ,2 bursts is likely to be dominated by contamination from lower redshifts. As such, 
the demographics (and indeed even the existence) of z > 8 bursts is largely unknown. To 
a lesser degree, the same problems arise for bursts from perhaps redshift 4 to 8. A typical 
consequence of this is that the GRB rate density cannot be rising as steeply out to very 
high redshifts as concluded by Fenimore & Ramirez- Ruiz (2000), Schaefer, Deng, & Band 
(2001), Lloyd- Ronning, Fryer, & Ramirez- Ruiz (2002), and Yonetoku et al. (2004). Our new 
conclusion effectively rejects the evidence that points to GRBs demonstrating a steady rise 
of the Universe's star formation rate to very high redshift. 

Is there any hope for using the BATSE bursts plus the luminosity indicators to get good 
demographics to moderate or high redshifts? Our experience is that the reliance on any one 
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luminosity indicator is poor. Even when combining five luminosity indicators the resultant 
error in redshift is 26%, and this still makes for substantial smearing of the underlying redshift 
distributions. However, at the 26% level, good work can be done with reasonable accuracy. 
Certainly, a detailed model of the smearing caused by the redshift uncertainties must be 
performed. Another requirement is that large number of BATSE bursts must be analyzed 
so as to avoid selection effects like in this paper. In all, we expect that a comprehensive 
study of BATSE bursts can pull out high redshift demographics provided the study (a) uses 
many of the luminosity indicators, (b) models the effects of the usual uncertainties, and (c) 
includes a large number of bursts without selection by the indicators. 

We thank Robert Mallozzi, now deceased, for his measures of E pea k, a, and j3. We 
acknowledge NASA grant NNG06GH07G for support of our work. 
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Table 1. Burst Properties and Bolometric Peak Flux 
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a This burst was chosen for its very high reported redshift (Band et al. 2004) based on the lag-luminosity relation. 

b This burst was chosen for its very high reported redshift (Fenimore & Ramirez-Ruiz 2000) based on the variability- 
luminosity relation. 

c This burst was chosen for its very high reported redshift (Yonetoku ct al. 2004) based on the E pea fe-luminosity relation. 
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Table 2. Burst Luminosity Indicators 
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Table 3. Burst Luminosity Relations 
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Table 4. Derived Redshifts and the One-Sigma Confidence Intervals 
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14.5 

1.6 

3.5 

3.0 

3.4 

2.0 

6.8 

5.8 

3.8 



2.3-20) 
1.9 - 18.5 ) 



1.4 - 
1.8 - 
( 3.4 

1.2 - 
1.4 - 
( 5.0 
2.1 - 
1.6 - 
1.8 - 
2.4 - 

1.3 - 
( 4.5 
1.8 - 

1.4 - 
0.5 - 
0.8 - 
5.0 - 
3.0 - 
1.6 - 

5.3 - 
1.6 - 
1.8 - 

( 3.2 
1.3 - 

1.5 - 
( 3.0 
0.9 - 

1.6 - 

1.0 - 

1.6 - 

1.1 - 

2.2 - 
2.2 - 

1.7 - 



6.9 ) 
18.7 ) 

- 20 ) 
7.4 ) 
7.5) 

- 20 ) 
20 ) 

10.7 ) 
18.2 ) 
20 ) 

8.6 ) 

- 20 ) 

18.5 ) 

7.7 ) 

I. 7) 

2.8 ) 
20 ) 
20 ) 
12.2 ) 

20 ) 
14.9 ) 
20 ) 

- 20 ) 
7.3 ) 

11.6 ) 

- 20 ) 

3.9 ) 

11.8 ) 
20 ) 

II. 6 ) 
5.0 ) 
20 ) 
20 ) 
19.7 ) 



0.9 (0.5-1.6) 

2.4 ( 1.3 - 5.4 ) 

3.1 ( 1.6 - 
0.6 ( 0.3 - 

11.6 ( 2.7 

2.2 ( 1.1 - 

16.7 ( 6.2 
3.1 ( 1.5 - 
20 ( 11.7 
11.5 ( 4.5 
1.6 ( 0.9 - 

9.5 ( 3.8 - 
5.8 ( 2.5 - 

8.3 ( 2.7 - 
0.9 ( 0.5 - 

16.8 ( 6.1 
7.5 ( 3.2 - 

4.4 ( 2.0 - 
5.4 ( 2.3 - 
19.4 ( 6.1 
5.8 ( 2.5 - 
20 ( 5.0 - 



6.8 ) 

1.1 ) 

- 20 ) 

4.9 ) 

- 20 ) 
9.6 ) 

- 20 ) 

- 20 ) 

3.2 ) 
20 ) 
18.5 ) 
20 ) 

I. 6) 

- 20 ) 
20 ) 

II. 2 ) 
20 ) 

- 20 ) 
17.7 ) 

20 ) 



0.7 ( 0.4 - 

10.3 ( 3.8 - 
17.6 ( 3.2 ■ 

1.2 ( 0.6 ■ 
2.9 ( 1.4 - 
3.5 ( 1.5 ■ 
0.6 ( 0.4 - 

4.3 ( 2.1 - 
20 ( 6.2 ■ 

14.4 ( 4.9 - 20 ) 
6.8 ( 2.7 - 20 ) 
3.2 ( 1.4 - 10.6 ) 
0.8 ( 0.5 - 1.5 ) 

3.4 ( 1.6 - 8.7 ) 



1.2) 
20 ) 
20 ) 
2.4 ) 
6.9 ) 
16.1 ) 
1.1 ) 
10.3 ) 
20 ) 



> 0.0 

> 0.0 

> 0.0 

> 0.0 

> 1.5 

> 0.5 

> 1.8 

> 0.5 

> 1.2 

> 1.4 

> 1.2 

> 0.6 

> 1.3 

> 3.1 

> 0.6 

> 1.3 

> 0.0 

> 0.8 

> 1.0 

> 3.0 

> 1.6 

> 1.5 

> 0.0 

> 0.9 

> 1.4 

> 0.0 

> 0.8 

> 0.3 

> 0.0 

> 0.9 

> 0.5 

> 1.6 

> 0.9 

> 0.2 

> 0.7 

> 1.6 



1.4 
3.7 
3.0 
1.8 
8.5 
2.4 



1.0-2 
2.3 - 
1.9 - 
1.2 - 
5.2 - 
1.5 - 



14.5 ( 8.2 

6.4 ( 3.4 - 
20 ( 14.9 
8.2 ( 5.1 - 
3.1 ( 2.2 - 
8.9 ( 4.5 - 

2.6 ( 1.8 - 

13.8 ( 5.4 

1.5 ( 1.1 - 
11.5 ( 5.6 
2.8 ( 1.7 - 

3.7 ( 2.3 - 

11.9 ( 5.4 
9.1 ( 5.6 - 
11.1 ( 5.6 
20 ( 8.9 - 
1.1 ( 0.8 - 

7.6 ( 4.8 - 



6.4 
1.3 
2.8 
4.7 
1.2 
3.9 



3.4 
1.0 
1.7 
2.7 
0.9 
2.4 



20 ( 8.3 - 
13.3 ( 6.0 
5.1 ( 2.8 - 
2.5 ( 1.7 - 
2.9 ( 1.9 - 
4.7 ( 3.1 - 



0) 

6.1) 
5.0 ) 
2.7 ) 
20 ) 
4.0 ) 

- 20 ) 
19.1 ) 

- 20 ) 
14.0 ) 
4.6 ) 
20 ) 
3.8) 

- 20 ) 

2.0 ) 

- 20 ) 
4.9 ) 
6.9 ) 

- 20 ) 
16.4 ) 

- 20 ) 
20 ) 

1.5) 
14.6 ) 
15.4 ) 
1.8) 
4.9 ) 
9.6 ) 

I. 7) 

7.1 ) 
20 ) 

- 20 ) 

II. 3 ) 
3.9 ) 
4.5 ) 
7.9 ) 



1.0 

2.3 

1.9 

1.2 

5.2 

1.5 

8.2 

3.4 

14.9 

5.1 

2.2 

4.5 

1.8 

5.4 

1.1 

5.6 

1.7 

2.3 

5.4 

5.6 

5.6 

8.9 

0.8 

4.8 

3.4 

1.0 

1.7 

2.7 

0.9 

2.4 

8.3 

6.0 

2.8 

1.7 

1.9 

3.1 



